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Abstract. A generalized version of the Spherical Infall Model (SIM) is used to study the effect of angular momentum 
on the final density profile of a spherical structure. The numerical method presented is able to handle a variety 
of initial density profiles (scale or not scale free) and no assumption of self-similar evolution is required. The 
realistic initial overdensity profiles used are derived by a CDM power spectrum. We show that the amount of 
angular momentum and the initial overdensity profile affect the slope of the final density profile at the inner 
regions. Thus, a larger amount of angular momentum or shallower initial overdensity profiles lead to shallower 
final density profiles at the inner regions. On the other hand, the slope at the outer regions is not affected by the 
amount of angular momentum and has an almost constant value equal to that predicted in the radial collapse 
case. 
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1. Introduction 

It is likely that dark matter halos are formed by the evolu- 
tion of small density perturbations in the early Universe. 
The matter contained in a perturbed region progressively 
detaches from the general flow and after reaching a ra- 
dius of maximum expansion it collapses to form an in- 
dividual structure. The most simple case is when this re- 
gion is spherical, isolated and undergoing a radial collapse. 
This is the spherical infall model (hereafter SIM). SIM has 
been extensively discussed in the literature (Gunn & Gott 
1972] ; Gott p7§ G unn pTT] ; Fillmore fc G oldreich |l984 



Frenk et al. |l988| ; Dubinski & Galberg |l99l[ Cr one et al 



1994; Navarro et al. 1997; Cole fc Lacey 1996; Husset 



al. 199E; Fukushinge & Makino 1997; Moore et al. 1998 



Jing & Suto 2000 ) showed that the profile of relaxed ha- 
los steepens monotonically with radius. The logarithmic 
slope a — — 4 1" p is less than 2 near the center and larger 
than 2 near the virial radius of the system. The value of 
a near the center of the halo is not yet known. Navarro 
et al. ( |l997j ) claimed a — 1 while Kravtsov et al. fll998| ) 
initially claimed a ~ 0.7 but in their revised conclusions 



Bertschinger 1985; Hoffm an & Shaham 1985, hereafter HS; 
White & Zaritsky |l992|) . 

The final density profile after a collisionless evolution 
of the matter depends on its initial density profile and the 
underlying cosmology. Self-similarity solutions (Fillmore 
& Goldreich |1984| Bertschinger |1985|, HS) show that a 



power-law initial density profile relaxes to a final density 
profile given by p{r) oc r~ a with a > 2. Furthermore, 
recent numerical studies that relax the assumption of self- 
similarity, also give final density profiles steeper than r~ 2 . 
Lokas & Hoffman ( 2000b ) found values of a in the range 
2 to 2.3. 

The density profiles of galactic halos do not se em to 
follow power laws. Numerical studies (Quinn et al. 1986; 



(Klypin et al. 200C ) they argue that the inner slope varies 
from 1 to 1.5. Moore et al. (1998) found a slope a = 1.5 
at the inner regions of their N-body systems. 

In this paper we study the final density profiles pre- 
dicted by the SIM, when non-radial motions are included. 
Studies concerning the role of non-radial motions have 
been presented by Ryden & Gunn (1987), Ryden (1988), 
Gurcv ich & Zybin ( |1988a|, |1988b| ), Avila-Reese et al. 
( |1998D , White & Zanteky (|1992j), Sikivie et al. Ql997j ) and 



recently by Nusser ( [200l|) . 

In Sect. 2 we discuss the SIM and the associated prob- 
lems. In Sect. 3 the description of the numerical method 
as well as the way the angular momentum is included is 
given. The initial conditions and the results are given in 
Sect. 4 and are summarized in Sect. 5. 
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2. The spherical infall model (SIM) 

SIM i s based on the physical pro cess d escribed in Gunn 
( 1977 ) and in Zaroubi & Hoffman ( 1993 ): In an expanding 
spherical region the maximum expansion radius £ (apap- 
sis) of a shell is a monotonic increasing function of its 
initial radius Xi and is given by the relation: 



C = 9{%i) = 



1 - fir 1 + Ai( Xi ) 



(i) 



where fi^ is the initial value of the density parameter of 
the Universe and Aj is the relative excess of mass inside 
the sphere of radius Xi, given by 



Ai(xi) 



M(xj) - M b { Xl ) 
M b ( Xi ) 



Xi JQ 



x 2 Si(x)dx. 



(2) 



In (|2|), M is the mass of the spherical region, M 
is the mass of the unperturbed Universe and Si is the 
spherically symmetric perturbation of the density field 
(Si(x) = P ^p~ Pt '' where p is the density and p oi is the 



constant density of the homogeneous Universe at the ini- 
tial conditions). The time spent for a shell to reach its 
above turnaround radius is: 



tta 



1 + Ai(xi) 



-IT, 



(3) 



where Hi is the value of Hubble's constant at the epoch 
of the initial conditions. The above equations are valid 
for bound shells, so the condition 1 — fi^ 1 + Aj(xi) > 
is satisfied. A shell, after reaching its turnaround radius, 
collapses, re-expands to a new (smaller) turnaround ra- 
dius and so on. The limiting value of this radius, after 
a large number of such oscillations, is the final radius of 
the shell, corresponding to the relaxed state of the sys- 
tem. In any reasonable potential a shell spends most of 



its orbital time near the maximum radius (Gunn 1977) 



So, if shell crossing -occurring during the collapse stage- 
had no dynamical consequences, the final distribution of 
mass could be approximated by the distribution resulting 
if every shell stopped at its maximum expansion radius. 
This should lead to a "turnaround density profile" pt a of 
the form 



<?) ''' :: - r ' ; (d^ 



(4) 



In deriving (Q) the conservation of mass is used (M(xi) — 
M(()). This is an important relation since this distribu- 
tion of mass is used as the initial one in SIM. SIM assumes 
that the collapse is gentle enough. This means that the or- 
bital period of the inner shell is much smaller than the col- 
lapse time of the outer shells (Zaroubi & Hoffman 1993[ ). 
This implies that the radial action <fv(r)dr, where v is 
the radial velocity, is an adiabatic invariant of the inner 
shell. As the outer shells collapse, the potential changes 
slowly and because of the above adiabatic invariant, the 
inner shell shrinks. The collapse factor depends on the 



time the mass of the outer shells (passing momentarily) 
spends inside the inner shell. Consider a shell with apap- 
sis C and initial radius Xi. The mass inside radius £ is a 
sum of two components. The first one, (permanent com- 
ponent, M p ), is due to the shells with apapsis smaller than 
C and the second (additional mass, M ad d) is the contribu- 
tion of the outer shells passing momentarily through the 
shell £. Because of the mass conservation, the permanent 
component is given by the following relation 



Mp(C) = M{xi) = -np bti xl[l + Ai{xi)]. 
The additional component is: 



M add (C) 



C 



R , dM(x) , 



(5) 



(6) 



In (g) , R stands for the radius of the system (the apapsis of 
the outer shell) and the distribution of mass M[x) is given 
by (0). P^(x) is the probability of finding the shell with 
apapsis x inside radius £, calculated as the ratio of the 
time the outer shell (with apapsis x) spends inside radius 
Q to its period. In the general case of non-radial collapse 
this ratio is given by the relation P((x) = I(()/I(x) where 



I(r) = 



dn 

v x (n) ' 



(7) 



where x p is the pericenter of the shell with apapsis x and 
v x {n) is the radial velocity of the shell with apapsis x as it 
passes from radius n. If the collapse is radial then x p = 0. 
After the calculation of M add the collapse factor f(xi) of 
a shell with initial radius Xi and apapsis £ — g(xi) is given 

by 



M P (C) 



M P (C) + M add {Q' 



(8) 



The final radius of the shell is x = f(xi)£ and mass con- 
servation leads to the following final density profile 



Pi?) =Aa(C)/ fri) 



din f(xi 



din g{xi) 



(9) 



The potential energy of the system Wf in the relaxed 
state is related to its total energy by the virial theorem 
Wf = IE. The energy of the system in a radial collapse 
is that at the turnaround epoch when all shells are as- 
sumed to have zero velocities at the same time. Thus 
Wf = 2W ta - A simple collapse factor satisfying this re- 
quirement is / = 0.5 for every shell, leading to similarity 
solutions (a final density profile parallel in log-log space to 
that of the turnaround epoch). However, the collapse fac- 
tor is not constant. N-body simulations (e.g. Voglis et al. 



1995, hereafter VHH) show that / is an increasing func- 



tion of the initial radius (or the turnaround radius) of the 
shell and its form is related to the initial profile of the 
density perturbation. 

In a radial collapse case eqn.(||) gives / — » as Xi — > 
resulting in very condensed central regions with very steep 
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density profiles (Lokas 2000a| ). As it is shown by Lokas 
& Hoffman (2000b) the inner slope of the density profile 
is between 2 and 2.3 even for low initial density peaks, 
far from the values derived from N-body simulations that 
range from 1 to 1.5 for the inner regions. On the other 
hand, the collapse factor in N-body simulations is hardly 
less than 0.05, even in the very central regions (VHH). 
This difference is easy to understand because no radial col- 
lapse exists in N-body simulations. During the expansion 
and the early stage of collapse particles acquire random 
velocities that prevent them from penetrating the inner 
regions. The consequence is the reduction of M a dd in (§) 
that leads to larger values of /. In this sense the approxi- 
mations based on constant / may be closer to the results 
of N-body simulations. The fit of the NFW profile is char- 
acteristic using an approximation given by del Popolo et 
al. ( |2000D . However, in such an approximation the final 
density profile depends only on the "turnaround density 
profile" and not on the amount of angular momentum ac- 
quired by the system during its expansion phase. 

3. The numerical method and the angular 
momentum 

The calculation of the collapse factor requires the evalu- 
ation of the integral in (^|) . Changing the variables from 
the turnaround radius to the initial one. this is written: 



M add (C) = 47rp M / P Xi {x'i)\ l + SM^xfdx't, 



where P Xi [x'^) = I(xi)/I(x' i ) with 



J(r) 



1 



dg(n) 



v g(x>)(g(n)) 



dn 



m . 



(10) 



(11) 



£ = g{xi), x' p — g~ 1 (x p ) where x p is the pericenter of 
the shell with initial radius x[ . The upper limit x b of the 
integral in (|Io|), is taken to be the initial radius of the 
sphere that has collapsed at the present epoch. The radial 
velocity v of a shell with apapsis x = g(xi) as it reaches 
the radius r = g(ri) is given by the conservation of the 
energy of the shell and is: 



^( r ) = 2[*(r)- £ J 



(12) 



where ^ equals minus the potential <i>,s x equals minus the 
specific energy and j x is the specific angular momentum 
of the shell. The potential '5 after the change of variables 
is given by the expression: 



*\g(ri)] = 



GM(x b ) 

g(xb) 



G 



M{xj) dg{x t ) 
g 2 {xi) dxi 



dxi, 



(13) 



where the distribution of mass M(xi) is that at the initial 
conditions. The energy of the shell is calculated by: 



ng(xi)} 



Jx 



2g 2 (x t ) 



(14) 



The angular momentum is introduced by the following 
scheme: 

Each shell expands radially from its initial radius Xi 
up to its maximum expansion radius x. At this stage a 
specific angular momentum j x is added, given by j x = 
L\J M{x)x = C-J M(xi)g(xi), where C is a constant. This 
way of introducing angular momentum is consistent with 
the angular momentum distribution in N-body simula- 
tions (e.g, Barnes & Efstathiou 1987) and does not in- 
troduce any additional physical scale. It has been used by 



Avila-Reese et al. ( 1998 ) and recently by Nusser ( 2001 ). 

In this way the apocenter r a of a shell is its turnaround 
radius x while its pericenter r p is found by the solution of 
the equation: 



2r 2 {*(r)-e x ]-j 2 x = 0. 



(15) 



The change of variables described above requires the so- 
lution for r of the equation 

2g 2 (r)[^(g(r))-e x ]-£=0, (16) 
where the potential \& is that of the turnaround epoch. 



Nusser (2001) proved the following two important 
properties: 

1. If the angular momentum is introduced in the above 
described way in a spherical system with a power law 
density profile, then all shells have the same eccentric- 
ity. In fact, if if = r p /r a then the following equation 
holds; 



GC 2 (<p- 2 - 1) = 



2[V(r p ) - fr(r a )]r a 
M(r a ) 



(17) 



The right hand side of the above equation can be ex- 
pressed in terms of (p in the case of a power law density 
profile and completes the proof. 
2. If the potential evolves adiabatically, given at time t by 
the relation ^t(r) = k(t)^>(r) with k(t) a slowly vary- 
ing function of t and 'J the potential at the turnaround 
epoch, and the radial action is indeed invariant, then 
the eccentricity of every shell remains constant during 
the evolution. 

The radial action can be written in the form: 



i^ 2 1) 



V(ur a ) - 1>(r a ) 
*(r p ) - #(r a ) 



(1 



J r =jX 

du. (18) 



In the case of a power law density profile the quantity 
Wf~i~WFy i s written in terms of <p and u. Since J r 
is constant, then tp is constant. 



Nusser (2001) used these properties to estimate the 
asymptotic behavior of the density profile near and far 
from the center of a system with a power-law initial den- 
sity profile. Unfortunately these properties do not hold for 
more realistic density profiles. Similar power -law density 
profiles have been used by Sikivie et al. ( |1997| ) who used a 
CDM power spectrum to estimate an "effective" exponent 



4 



N. Hiotelis: Density profiles in a spherical infall model with non-radial motions 




1.5 2 
x (Mpc) 

Fig. 1. Overdensity profiles versus the distance x from a 
3tr peak. Dotted curve corresponds to a smoothing length 
of 1.2Mpc while the dashed and solid curves correspond 
to smoothing lengths of 0.6 and 0.3Mpc respectively. The 
values of < 8(x) > and x are both normalized to the 
present. 



for this power-law on the galactic scale. Sikivie et al. use 
a self-similar evolution of the system in order to calcu- 
late its properties at the relaxed state. Unfortunately, the 
self-similar evolution is not valid for more realistic initial 
density profiles (non power-law profiles). The numerical 
method presented in our study is more general. It is able to 
deal with various initial density profiles (scale or not scale 
free) and the assumption of self-similarity is not required. 
Our results (presented in Sect. 4) have derived for realistic 
initial density profiles that have a finite value of the den- 
sity perturbation at the location of the peak. Therefore, 
these results could be more reliable at least regarding the 
final state of the central region of the system. Moreover, it 
is shown that taking into account the angular momentum, 
final density profiles are well fitted by two power law den- 
sity profiles with slopes less than 2 at the central regions 
of the systems and larger than 2 at the outer regions. This 
class of final density profiles is consistent with the results 
of N-body simulations (eg Subramanian et al. ( pOOO )). 



4. The initial density profile and the results 

The averaged overdensity profile < S(x) > at distance x 
from a na = n^ 1 / 2 (0) extremum of a smoothed density 
is given in Bardeen et al. (1986, hereafter BBKS) by the 
equation: 



:[7 ?W 



7(1-78)' 



< S(x) > 
V 2 £(z)]/£(0) 1/2 , 



n£(x) 



(19) 



where 7 = 7(4)/[/(2)/(6)] 1 / 2 and R* = [3/(4)//(6)] 1 / 2 
with 1(1) = J °° k l P(k)dk. 

In the above relations £ is the correlation function and 
P the power spectrum. These are related by: 



2tt 2 x 



P(k)ksin(kx)dk 



The function #(717, 7) is given by the relation 



0(n7,7) 



T\ 2 > 



3(1 - 7 2 ) + (1.216 - 0.9 7 4 )e- 
[3(1 - 7 2)+ 0.45 + (^) 2 ] I + ^ 



AM 



in the range 1 < n 7 < 3. The rms mass excess, -p-, 
a sphere of radius x, is given by: 



1 



(2tt 2 )^ x 6 



P(fc)(sin kx — kx cos kx)' 



dk 



(20) 

(21) 
within 

(22) 



We used the spectrum calculated by BBKS for a CDM- 
dominated Universe with fl — 1 and h = 0.5. This is given 
by the equation 



Pcdm = Afc-^Ml + 4.164fc)] 2 [G(fc)]^ 

where G{k) is the following polynomial of k. 

GQs) = 192.9 + 1340fc + 1.599 x 10 5 fc 2 
+1.78 x 10 5 /c 3 + 3.995 x 10 6 fc 4 



(23) 



(24) 



The above spectrum is smoothed on various scales accord- 
ing to the relation: 



P(k) = PcDAie 



(25) 



Regarding the smoothing of the above spectrum three 
cases are examined. In case A the spectrum is smoothed on 



cl SCctlc kf 



1.2 Mpc, in case B on a scale k c 



0.6 



Mpc and in case C on k c 1 = 0.3 Mpc. The first scale 
length corresponds to a mass 0.5 x 10 12 Mq the second to 
a mass 6.25 x 1O 1O M and the third to 7.8 x 10 9 Af Q . The 
constant of proportionality A is chosen by the condition 
(Ts = 0.7, (the rms mass excess within 8 Mpc to be 0.7). 
Then < 6(x) > is calculated for n = 3 and plotted in 
Fig. |l|. The dotted curve corresponds to the case A, the 
dashed curve to the case B, while the solid one to the case 
C. We note that the values presented in this figure are 
normalized to the present. 

In the case of linear growth of the overdensity the fol- 
lowing hold (Gunn & Gott 1972): a shell with initial ve- 
locity equal to the Hubble flow and an initial comoving 
radius x has expanded up to a maximum radius r max in a 
time tta (given by (^|), fli = 1 in our case). This maximum 
radius is given by the relation 



x/A, 



(26) 



while the collapse time of the shell, t c , is related to the age 
of the Universe,^, by 



(27) 
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However calculating 

3 r 

A = — / < 6(u) > u 2 du, 
x Jo 



(28) 



and using ( |27| ) and ([28|) the collapse time of a shell and its 
turnaround radius are found. In fact the condition t c /to = 
1 gives the value of A for the shell that collapses today. 
Then x is found by solving numerically ( p8| ) and r max is 



calculated by (26). In case A the mass inside the shell that 
collapses today is about 9.8 x 10 12 Mq and the radius of 
maximum expansion is 1130 Kpc. For the case B the mass 
is 4.3 x 10 12 Mq and the radius of maximum expansion is 
870 Kpc while for the case C the values are 2.1 x 10 12 Mq 
and 685 Kpc respectively. The values of x resulting from 
the numerical solution of ( |28| ) are 3.18, 2.45 and 1.93 Mpc 
respectively. Finally, the initial conditions at redshift zi 
can be derived by dividing both x and < S(x) > by 1 + Zj. 
Using 1 + Zi= 1000 the value of x b in @ are 3.18, 2.45 
and 1.93 Kpc for the cases A, B and C respectively. 

The amount of the angular momentum in the system 
is adjusted by the value of C, (see Sect. 3), and is measured 
by the value of the dimensionlcss spin parameter 



A 



L\E\i 



(29) 



where L, E and M are the total angular momentum, the 
total energy and the total mass of the system respectively. 
The mean value of A resulting from N-body simulations 
(e.g. Efstathiou & Jones |1979| , Barnes & Efstathiou [1987D 
seems to be about 0.05. In our calculations we used differ- 
ent values of A with a maximum of 0.12. The maximum 
value of A corresponds to C — 0.26. 

The resulting density profiles are fitted by a two-power 
law curve of the form 



Pfit(r) 



Pc 



where the fitting parameters p c , r s , p and \x are calculated 
finding the minimum of the sum 



(\0g ns i ,,(/',: ) 

1=1 



log Pfitin) 



log psim{ti) 



(31) 



where psiM are the predictions of SIM. NP is the num- 
ber of points where the density is found. We used 100 
points equally spaced on a log scale. The estimation of the 
above fitting parameters is done using the unconstrained 
minimizing subroutine ZXMWD of IMSL mathematical 
library. The quality of fit is very good as can be seen in 
the following three Figs. 

In Fig. U the final density profiles for the case A are 
shown for three different values of A. The radial collapse 
(A = 0) corresponds to the higher solid curve. The inter- 
mediate curve corresponds to A = 0.05 while the lower 
one to A = 0.09. Dotted lines are the fits by the above de- 
scribed two-power law density profile. Distances are nor- 
malized to the virial radius r vir which is 448 kpc while 
densities are normalized to the critical density p cr a- The 




-2 -1.5 -1 
log(r/r vir ) 

Fig. 2. Case A. Final density profiles derived for three 
different values of the spin parameter A. Solid curves: SIM 
predictions. Dotted curves: the fits by a two-power law. 
From the bottom of the figure the curves correspond to 
A = 0.09, 0.05 and 0.0 (radial collapse) respectively. 



(30) - 




-3 



-2.5 



-2 -1.5 -1 

log(r/r.) 



-0.5 



Fig. 3. Case B. As in Fig.^. Solid curves: Density profiles 
derived from the SIM. From the higher to the lower curve 
the values of the spin parameter are 0.0, 0.05, 0.09 and 0.12 
respectively. Dotted curves: The fits of the solid curves by 
a two-power law density profile. 



virial radius is the radius of the sphere with mean density 
« 178 ti mes th e present density of the Universe p cr it (Cole 
& Lacey 1996 ). The virial mass of the system (the mass 
contained inside the virial radius) is about 4.7 x 10 12 Mq. 
It is characteristic that the efficiency of the angular mo- 
mentum leads to shallower density profiles in the inner 
regions of the system. Additionally at the outer regions 
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-3 -2.5 -2 -1.5 -1 -0.5 

log ( r / r . ) 

Fig. 4. Case C. As in Fig.||. From the higher to the lower 
solid line the values of A are 0.0,0.05,0.09 and 0.12 re- 
spectively. Dotted lines are the fits of the solid lines by a 
two-power law density profile. 

the density profile does not change even for the maxi- 
mum amount of angular momentum used. We note that 
for larger values of the spin parameters the density profile 
becomes unrealistic (increasing at the inner regions). As 
it will be shown below this is a consequence of the shallow 
initial profile of this case. 

The results for the case B are presented in Fig. [| The 
values of the spin parameter are 0., 0.05, 0.09, and 0.12. 
The virial radius of the system is about 347 Kpcs and 
contains a mass of about 2.16 x 10 12 Mq. 

The case C is presented in Fig. |] for the same values 
of A as in the case B. The virial radius of this system is 
273 Kpc and its virial mass 1.06 x 10 12 Mq. 

Note that in the cases B and C the profile of the density 
decreases even for larger values of A than used in case A, 
because of their steeper initial density profiles. 

In the following three figures the collapse factors for 
each case are presented. Fig. || shows the collapse factor 
/ of mass M on a logarithmic scale. The role of angular 
momentum is clear. Larger values of A lead to smaller 
values of / and consequently to shallower density profiles 
as shown in Figs 0, || and [|. Fig.^ refers to case B and Fig. [7] 
to case C. It is clearly shown in the above three figures 
that the collapse factor at the outer regions of the system 
is not affected by the amount of the angular momentum 
and it is almost the same as that of the radial collapse case. 
The efficiency of angular momentum in creating shallow 
density profiles depends on the initial density profile. This 
can be shown in the following three Figs where the slope 
of the two-power law fit versus radius is plotted. This is 
given by the relation 

1"s 




_q ill lllll llllll llllll ill lllll llllll 11 

-3 -2.5 -2 -1.5 -1 -0.5 0.5 1 
logM(10 12 M o ) 

Fig. 5. The collapse factor of mass M versus M on a 
logarithmic scale for case A. From the lower line the values 
of A are 0.0, 0.05, 0.09 respectively. 

Fig. H corresponds to case A. The higher line is the slope 
resulting after a radial collapse (A = 0.0) where the val- 
ues of a in the interval 0.001r^j r to the range 
2. to 2.25. The intermediate line corresponds to A = 0.05 
while the lower line corresponds to A = 0.09. Fig.|| shows 
the slope for the case B. The lines, from the higher to the 
lower, correspond to A = 0.0,0.05,0.09 and 0.12 respec- 
tively. The results of case C are shown in Fig.[l0| for the 
same values of A as in case B. The values of a for the three 
cases at r = 0.0017W are clearly shown in the above fig- 
ures. At r = 0.01r vir and for A = 0.05 the values of a are 
1.76, 1.82 and 1.83 for the cases A, B and C respectively. 
A similar trend for the inner regions -smaller a for shal- 
lower initial density profile-is also clear for all values of 
A. This trend is reversed at the outer regions. The slopes 
at r = r V i r are approximately 2.25, 2.20 and 2.15 for the 
three cases respectively. It is also clear from Fig.[l(] that 
the radial collapse case leads to an almost exact power law 
profile. In this case the slope at r = 0.01r„i r is 2.12 while 

at T — — Tyir is 2.15. 

5. Conclusions 

The predictions of the SIM presented in this paper are 
summarized as follows: 

1. Radial collapse does not lead to power law final den- 
sity profiles. However the difference in the slopes be- 
tween the inner and the outer regions of the system is 
not large. Decreasing the smoothing scale of the power 
spectrum (leading to steeper initial profiles) this differ- 
ence becomes smaller, leading to an almost power law 
for steep enough initial density profiles. The slopes in 
the radial collapse case are, in agreement with theoret- 
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Fig. 6. As in Fig.| but for case B. The values of A, from Fig. 8. The slope of the final density profile versus ra- 
the lower to the higher line are 0.0,0.05,0.09, and 0.12 dius for case A. Distance is normalized to the virial ra- 
dius. From the lower to the higher line the values of A are 
0.09, 0.05 and 0.0 respectively. 




log M (10 ,z M o ) 
Fig. 7. As in FigJ| but for the case C. 




ical predictions (HS, Bertschinger 1985| ), in the range 
2 to 2.25. 

Angular momentum leads to shallower inner density 
profiles. The inner slope depends on the amount of the 
angular momentum, measured in our results by the 
value of the spin parameter, and on the form of the ini- 
tial density profile. Angular momentum becomes more 
efficient, in decreasing a, for shallower initial density 
profiles. 

The slope of density profiles does not change signifi- 
cantly at the outer regions of the system even in cases 
where a large amount of angular momentum is as- 



log(r/r vir ) 

Fig. 9. As in FigJ| but for the case B. From the lower to 
the higher line the values of A are 0.12,0.09,0.05 and 0.0 
respectively. 



signed to the system. At r = r V i r the slope is approxi- 
mately that of the radial collapse case. 

We note that the above results are limited by a large 
number of assumptions, by the specific underlying cosmol- 
ogy and the particular form of the power spectrum used. 
However, they show systematic trends that could help us 
to better understand the relation between the initial con- 
ditions and the final density profiles. If things go the way 
described above, then the results of N-body simulations 
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could be approximated by adding angular momentum to 
a case where the radial collapse results in a r~ 3 density 
profile at the outer regions. However, the role of different 
parameters of the problem is under study. 
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